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Abstract. An analysis of the equations used for modeling thermal arsenic diffusion in silicon has been 
carried out. It was shown that for arsenic diffusion governed by the vacancy-impurity pairs and the 
pairs formed due to interaction of impurity atoms with silicon self-interstitials in a neutral charge state, 
the doping process can be described by the Fick’s second law equation with a single effective diffusion 
coefficient which takes into account two impurity flows arising due to interaction of arsenic atoms with 
vacancies and silicon self-interstitials, respectively. 

Arsenic concentration profiles calculated with the use of the effective diffusivity agree well with ex¬ 
perimental data if the maximal impurity concentration is near the intrinsic carrier concentration. On 
the other hand, for higher impurity concentrations a certain deviation in the local regions of arsenic 
distribution is observed. The difference from the experiment can occur due to the incorrect use of effec¬ 
tive diffusivity for the description of two different impurity flows or due to the formation of nonuniform 
distributions of neutral vacancies and neutral self-interstitials in heavily doped silicon layers. 

1 Introduction 

Arsenic is the basic n-type impurity in the technology of manufacturing silicon semiconductor devices 
and integrated microcircuits. It is due to the small difference between the atomic radii of arsenic and 
silicon that results in a high value of arsenic solubility limit in silicon. At present, arsenic doped regions 
are formed by ion implantation with subsequent thermal annealing of implanted layers. The thermal 
treatment provides the electric activation of impurity atoms, recrystallization of amorphous layers, and 
annealing of radiation defects created by ion implantation. During annealing a transient enhanced diffu¬ 
sion of impurity atoms occurs that results in the increase in the doped region dimensions and, respectively, 
in the increase of the p-n junction depth. On the other hand, due to the strong dependence of the arsenic 
diffusivity on impurity concentration, there occurs the formation of the specific impurity distribution of a 
“box-like” form. These impurity distributions are characterized by a practically constant arsenic concen¬ 
tration within the limits of the doped region and a very sharp fall of impurity concentration in the bulk 
of a semiconductor. Such kind of impurity atom distribution is very important in semiconductor manu¬ 
facturing because it imparts high electrophysical characteristics to semiconductor devices. It follows from 
the presented analysis that knowing the concentration dependence of arsenic diffusivity is very important 
for solving the problems of doping process modeling in the semiconductor technology. Investigation of 
different concentration dependences of arsenic diffusivity is the goal of this research. 
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2 Model of diffusion of arsenic atoms 


To investigate the concentration dependence of arsenic diffusivity, it is reasonable to consider the diffusion 
processes which are independent of the factors that complicate the diffusion description, such as arsenic 
clustering or precipitation, annealing of the defects created by ion implantation, oxidation of the surface 
of a semiconductor, etc. Therefore, let us consider thermal arsenic diffusion from a constant dopant 
source on a surface that ensures the maximal impurity concentration below a solubility limit of arsenic 
in silicon. Let us also assume that the distributions of neutral point defects responsible for impurity 
diffusion are approximately uniform within the diffusion zone il] and the local charge neutrality is valid 
mmm- Then, the system of equations for impurity diffusion due to the formation, migration, and 
dissociation of the “impurity atom — vacancy” and “impurity atom — silicon self-interstitial” pairs, 
being in a local thermodynamic equilibrium with the substitutionally dissolved impurity atoms and point 
defects participating in the pair formation [B] , takes the following form: 
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Here C and C B are the concentrations of substitutionally dissolved arsenic atoms and background 
impurity of the opposite type of conductivity, respectively; x is the concentration of electrons n, normal¬ 
ized to the intrinsic carrier concentration m\ h (C, C B ) is the factor describing the influence of built-in 
electric field on drift of charged pairs; D Fx , D F ~, and Df 2 ~ are the partial intrinsic diffusivities of 
dopant atoms due to interaction with neutral, singly, and doubly charged vacancies, respectively; /3 F 
and f3 F are the empirical constants that describe the relative contribution of singly and doubly charged 
vacancies to the arsenic diffusion. Likewise, D Fx , D F ~, and D F2 ~ are the partial intrinsic diffusivities of 
arsenic atoms due to interaction with neutral, singly, and doubly charged self-interstitials, respectively; 
/3 F and f3 F are the empirical constants that describe the relative contribution of singly and doubly charged 
self-interstitials to the impurity diffusion. 
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It follows from Eq. © and expressions ED, Q, ©, and ED that to calculate the arsenic concen¬ 
tration distribution it is necessary to know the intrinsic diffusivities Df and Df that take into account 
the interaction of impurity atoms with vacancies and self-interstitials, respectively, and 4 values of the 
empirical constants fif, fif, /3f, ff that describe the relative contribution of the charged point defects. 
Unfortunately, it is very difficult to find these values because even the values of the relative contribution 
of the diffusion mechanism via self-interstitial to the total arsenic transport 
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are characterized by well-marked differences as follows from [7j. For example, f 1 = 0.09 in the temperature 
range 975 1200 °C according to the data of (8], whereas f 1 = 0.35-0.55 for a temperature of 1000 °C 
according to [9]. 

In view of the wide scatter of the values of f 1 and the complexity of separate determination of the 
parameters /df, fif, /3f, and /3f, a more simple diffusion equation: 


— = V{ [D(x)] h (C,C B ) VU} , 

with the effective diffusivity 


( 8 ) 
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is very often used for modeling arsenic diffusion. The empirical constants and /?2 in expression ED 
are found from the best fit to experimental data. For example, in Ref. |10j a fitting routine based on 
different experimental data was used and temperature dependences of the arsenic partial diffusivities Df , 
Df , D~~ due to neutral, singly, and doubly charged intrinsic point defects were evaluated, which allows 
one to calculate the values of A = Df + Df + D^~ , fix = D~ / Df and /?2 = Df~/Df for different 
temperatures of thermal treatment. 

It is worth noting that impurity diffusion due to the “impurity atom silicon self-interstitial” pairs 
can occur basically via interaction with self-interstitials in the neutral charge state. Then, the description 
of impurity diffusion by Eq. ED and concentration dependence of diffusivity ED is identical to the use of 
the initial equation m and concentration dependence ED- Really, it is possible in this case to present 
the concentration dependence ED as a sum of concentration dependencies ED written for diffusion due to 
the vacancy-impurity pairs and for diffusion via the “impurity atom — silicon self-interstitial” pairs: 
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It follows from Eq. (1101) that 
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3 Simulation of arsenic diffusion in silicon 

In Fig. |T] the arsenic concentration profile calculated for the case of thermal diffusion from a constant 
source on the surface of a semiconductor is presented. The experimental data obtained in mi are used 
for comparison. In Ref. mi, arsenic diffusion was carried out at a temperature of 1108 °C for 120 min. 
The arsenic concentration profile was measured by neutron activation analysis. 

The analytical solution for the constant value of arsenic diffusivity, i.e., without taking account of the 
concentration dependence D = D (y) and the drift of charged pairs in built-in electric field, is 

C(x, t) = C s erfc (13) 

which is also presented by the dotted curve in Fig. |I] Here Cs = const is the arsenic concentration 
on the surface of a semiconductor, L = 2 \JDit is the characteristic length of impurity diffusion and 
erfc ( x ) = 1 — erf (x) is the additional error function. 



Depth (pm ) 

Figure 1: Arsenic concentration profiles formed by thermal diffusion from a constant source on the 
silicon surface. The dotted and solid curves represent the arsenic concentration profiles calculated from 
analytical solution m and obtained by numerical solution of the diffusion equation that takes account of 
the concentration dependence of arsenic diffusivity and drift of the charged species in the built-in electric 
field. Diffusion temperature is 1108 °C for 120 min. Filled circles are the experimental data of m 


The intrinsic diffusivity Di = 2.55xl0“ 6 pm 2 /s and the parameters /?i = 611 and fa = 44.94 that 
prescribe a concentration dependence of impurity diffusivity (J2|) at a temperature of 1108 °C have been 
obtained from the expression of [TO] : 
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It is clearly seen from Fig.[T]that the account of the drift of charged pairs (As + D x ) and (As + D 2_ ) and 
also the account of the concentration dependence of arsenic diffusivity due to the increased concentrations 
of charged defects D - and D 2 ~ in the region with higher impurity concentration results in a good 
agreement of the calculated profile with experimental data. It is worth noting that the analytical solution 
gives smaller concentration values and, respectively, smaller depth of the p-n junction in comparison with 
the experiment. 

Based on the good agreement of numerical calculations with experimental data m, one can carry 
out modeling of arsenic diffusion for higher impurity concentration at the surface of a semiconductor C$ 
= 2.0xl0 8 /im -3 . It follows from Ref. [T2] that the maximal value of equilibrium electron concentration 
n e that can be obtained in arsenic doped silicon at a temperature of 1108 °C is equal to 4.24xl0 8 ^m~ 3 . 
A further increase of electron concentration is limited by the clustering of impurity atoms. Taking into 
account these data on n e , it might be expected that for Cs = 2.0xlO 8 pm~ 3 the arsenic clustering is 
negligible and can be excluded in diffusion modeling. The arsenic concentration profile after thermal 
treatment for 30 min calculated under these assumptions is presented in Fig. [2] 



Depth (pm ) 

Figure 2: Arsenic concentration profiles formed by thermal diffusion from a constant source on the 
silicon surface. The dotted and solid curves represent the arsenic concentration profiles calculated from 
analytical solution m and obtained by numerical solution of the diffusion equation that takes account of 
the concentration dependence of arsenic diffusivity and drift of the charged species in the built-in electric 
field. Diffusion temperature is 1108 °C for 30 min 

It can be seen from Fig. [2] that the form of arsenic profile differs substantially from the distribution 
presented in Fig. |T|and from analytical solution (TT3l) . Indeed, due to the significant increase in arsenic 
diffusion in the high concentration region, the impurity distribution in this part of the profile becomes flat 
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and differs substantially from the erfc-function. On the other hand, the arsenic profile becomes abruptly 
falling in the bulk of a semiconductor and it is characterized by a significant concentration gradient. It 
means that using a strong nonlinear dependence of arsenic diffusivity on impurity concentration, it is 
possible to form concentration profiles of electrically active arsenic atoms with an ideal “box-like” form. 

Let us compare now the calculations of diffusion process characterized by a strong nonlinear concen¬ 
tration dependence of diffusivity with experimental data. Figure [3] presents the arsenic concentration 
profile calculated for the treatment temperature of 1050 °C and duration of 60 minutes. The parameters 
of arsenic diffusivity Di = 5.8134xl0 -7 /im 2 /s; /3\ = 238.9; and fa = 12.7 were taken from [TU]. The value 
of the impurity concentration at the surface of a semiconductor Cs has been chosen equal to 2.3xlO 8 
/im~ 3 . This value is below the maximal equilibrium electron concentration value n e = 3.566x10 s p m -3 
for the annealing temperature under consideration Therefore, it is possible to assume that arsenic 
clustering does not play an important role in diffusion. 
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Figure 3: Arsenic concentration profiles formed by thermal diffusion from a constant source on the 
silicon surface. The dotted and solid curves represent the arsenic concentration profiles calculated from 
analytical solution m and obtained by numerical solution of the diffusion equation that takes account of 
the concentration dependence of arsenic diffusivity and drift of the charged species in the built-in electric 
field. Diffusion temperature is 1050 °C for 60 min. Filled circles are the experimental data of m 

It can be seen from Fig. [3] that the calculated arsenic concentration profile agrees well with experimen¬ 
tal data both in the region of the high impurity concentration and in the region of the low concentration 
of impurity atoms including the region of p-n junction. On the other hand, some difference between the 
calculated and experimental profiles is observed in the region where a strong decrease in the impurity 
concentration begins. For comparison, in Fig. [4] the results of simulation of arsenic diffusion based on 
the model of m are presented. It is supposed in this model that doubly charged point defects do not 
significantly contribute to arsenic diffusion. The following parameters that describe arsenic diffusion have 
been used in simulation: Di = 5.537xl0 -7 pm 2 /s, /3i = 100, and /?2 =0 [13]. The impurity concentration 
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at the surface Cs has been chosen equal to 2.9x10 s /jm 



Figure 4: Arsenic concentration profiles formed by thermal diffusion from a constant source on the silicon 
surface. The dotted and solid curves represent the arsenic concentration profiles calculated from analytical 
solution m and obtained by numerical solution of the diffusion equation. Diffusion temperature is 1050 
°C for 60 min. Filled circles are the experimental data of m 


It can be seen from Fig. 0] that both the arsenic concentration profiles obtained by analytical solution 
and those by numerical calculation agree well with experimental data excluding the region of a low 
impurity concentration. It is worth noting that here the approximate analytical solution of the diffusion 
equation (O for /?2 =0 obtained in [Hj was used: 


where 


C(x, t ) 
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It was assumed in |14| that the concentration dependence of effective diffusivity is described by the 
following expression: 


DC {X ) h (C,C°) = 2 (A)\ 


(17) 
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4 Analysis of the results obtained by simulation of impurity 
diffusion with different concentration dependencies of arsenic 
diffusivity 

The arsenic concentration profile presented in Fig. |T] shows that the concentration dependence of arsenic 
diffusivity obtained in [TO] provides excellent agreement with experimental data on impurity concen¬ 
tration closed to the value of n;. On the other hand, if the arsenic concentration C » ni and the 
concentration profile of electrically active arsenic gets a “box-like” form (see Fig. [3]), a difference from the 
experimental data is observed in the local region where a strong decrease in the impurity concentration 
begins. Therefore, for solving the problems arising in modeling silicon doping with arsenic it seems more 
convenient to use the concentration dependence proposed in Il3l which is characterized by neglecting the 
interaction of impurity atoms with doubly charged point defects. Indeed, in this case, a good agreement 
with experiment is observed in the whole region of high impurity concentration and only a small difference 
takes place in the region of low arsenic concentration which is less important for subsequent calculations 
of electrophysical parameters. Taking into account the results of m, one can suppose that this difference 
appears due to the overlooked possibility of direct migration of arsenic interstitial atoms. 

Incomplete agreement of the calculated arsenic profile with the experimental ones can also be due to a 
certain inadequacy of using Eq. (HI) in a heavily doped region C » rii, instead of a more exact equation 
©■ Indeed, with increase in the impurity concentration, the concentration of electrons becomes greater 
than rij. As a result, the concentrations of negatively charged silicon interstitials are also increasing. 
If these defects play an important role in the impurity transport, the description of diffusion based 
on Eq. © becomes incomplete and it is necessary to use more exact equation ©. However, there 
can exist another explanation of the difference between numerical calculations and experimental data. 
Indeed, increase in the concentration of substitutionally dissolved impurity atoms as well as increase in 
the electron concentration can result in the formation of nonuniform distributions of neutral vacancies 
and silicon self-interstitials m • In this case, both Eq. © and Eq. © give inadequate description of 
the impurity diffusion. The solution of this problem requires further investigations. 

5 Conclusions 

The analysis of the equations used for modeling thermal arsenic diffusion in silicon has been carried out. 
With this purpose, different concentration dependencies of impurity diffusivity that are used for modeling 
arsenic transports processes have been analyzed. It is supposed that diffusion is carried out due to the 
formation, migration, and dissociation of “impurity atom — vacancy” and “impurity atom silicon self- 
interstitial” pairs which are being in a local thermodynamic equilibrium with substitutionally dissolved 
impurity atoms and point defects participating in pair formation. It is shown that for arsenic diffusion 
governed by the vacancy impurity pairs and by the pairs formed due to the interaction of impurity atoms 


with silicon self-interstitials in a neutral charge state, the doping process can be described by the Fick’s 
second law equation with a single effective diffusion coefficient which takes into account two impurity flows 
arising as a result of the interaction of arsenic atoms with vacancies and silicon interstitials, respectively. 

Arsenic concentration profiles calculated with the use of the effective diffusivity taking into account 
two different impurity flows agree well with experimental data if the maximal impurity concentration 
is close to the intrinsic carrier concentration n *. On the other hand, similar calculations for impurity 
concentrations above n, are characterized by a certain deviation from experimental data in local regions 
of arsenic distribution. So, the use of a diffusion model that takes account of the significant contribution 
of doubly negatively charged point defects leads to the deviation in the region where a significant decrease 
in the impurity concentration begins. On the other hand, the neglect of the contribution of these defects 
leads to a difference in the region of a low impurity concentration. This difference can be due to the 
incorrectness which arises if a single effective diffusivity is used for the description of two different impurity 
flows. On the other hand, the origin of these small differences can be the formation of nonuniform 
distributions of neutral vacancies and neutral silicon self-interstitials in heavily doped layers. 
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